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' Abstract 
, 

We start by formulating geometrically the Newton's law for a classical free par- 
ticle in terms of Riemannian geometry, as pattern for subsequent developments. In 
fact, we use this scheme for further generalisation devoted to a constrained particle, 
to a discrete system of several free and constrained particles. 

For constrained systems we have intrinsic and extrinsic viewpoints, with respect 
to the environmental space. In the second case, we obtain an explicit formula for 
the reaction force via the second fundamental form of the constrained configuration 
space. For multi-particle systems we describe geometrically the splitting related to 
the center of mass and relative velocities; in this way we emphasise the geometric 
source of classical formulas. 

Then, the above scheme is applied in detail to discrete rigid systems. We start 
by analysing the geometry of the rigid configuration space. In this way we recover 
the classical formula for the velocity of the rigid system via the parallelisation of 
Lie groups. Moreover, we study in detail the splitting of the tangent and cotangent 
environmental space into the three components of center of mass, of relative veloc- 
ities and of the orthogonal subspace. This splitting yields the classical components 
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of linear and angular momentum (which here arise from a purely geometric con- 
struction) and, moreover, a third non standard component. The third projection 
yields an explicit formula for the reaction force in the nodes of the rigid constraint. 
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Introduction 

The original approach to classical mechanics is based on the Newton's law. This is still 
used and popular mainly in the literature devoted to applied sciences and engineering, 
even if it is not very sophisticated from the mathematical and geometrical viewpoint (see, 
for instance, P US CSl 12^] ) • 

On the other hand, an approach to mechanics based on modern differential geom- 
etry has been developed and became more and more popular in the last decades. This 
viewpoint is achieved in terms of Riemannian, Lagrangian, Hamiltonian, symplectic, vari- 
ational and jet geometry. A very huge literature exists in this respect (see, for instance, 
PElEliainilSl[ini[IlllISlliniinilIHlliniiniE2ll21[24^ These methods have been very 
successful for understanding several theoretical aspects and for the solution of several 
concrete problems, and have stimulated a large number of further classical and quantum 
theories. 

In this paper, our aims are more specific and foundational. Namely, we reformulate 
classical mechanics of a system with a finite number of particles and rigid systems, in terms 
of the Newton's law, in a way which, on one hand, is closer to the classical treatment 
of the subjects and, on the other hand, is expressed through the modern language of 
differential geometry. 

Our approach is addressed both to differential geometers, who could easily get me- 
chanical concepts written in their language, and to mathematical physicists, who are 
interested in a mathematically rigourous foundation of mechanics. 

In fact, several ideas have been achieved independently by differential geometers and 
mathematical physicists in different contexts and with different purposes and languages. 
Sometimes, facts which appear in one of the two disciplines as easy and elementary may 
correspond to more difficult and fundamental facts in the other discipline. We believe 
that linking those facts provides a new insight on classical matters and yields new results 
as well. For the above reasons, from time to time, we recall some classical facts of one of 
the two areas which are possibly not very familiar to experts of the other area. 

Thus, this paper, in spite of the sophisticated mathematical language, in comparison 
to the standard literatures of mechanics, analyses concrete mechanical contents. 

On the other hand, this paper provides the classical background for a covariant ap- 
proach to the quantisation of a rigid body, which is the subject of a subsequent paper 

m 

The guideline of our approach is the description of mechanics of a system of n free and 
constrained particles, including a rigid system, in terms of the Riemannian formulation 
of mechanics of one particle. 

We start by recalling the mechanics of one free particle moving in an affine Euclidean 
configuration space. We express the Newton's law in terms of covariant derivative. In 
several respects, it is convenient to introduce forces as forms (instead as vector fields) 
from the very beginning. 

Then, we can naturally apply this Riemannian approach to the mechanics of a con- 
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strained particle. We have an intrinsic and an extrinsic viewpoint related to the embed- 
ding of the constrained configuration space into the environmental space. In particular, 
we use the Gauss' Theorem concerning the splitting of the Riemannian connection in 
order to get an explicit expression of the reaction force via the 2nd fundamental form of 
the constrained configuration space. 

Next, we describe the mechanics of a system of n free particles, as one free particle 
moving in a higher dimensional product configuration space. For this purpose, it is neces- 
sary to introduce a weighted metric (besides the standard product metric). Of course, in 
the case on n free particles, we have the additional projection on the single particle spaces. 
Furthermore, we have the splitting of the configuration space into the affine component 
of the center of mass and the vector component of relative distances. The 1st splitting 
can be used to achieve information on the single particles and is orthogonal with respect 
to both metrics. The 2nd splitting has a fundamental role and is orthogonal only with 
respect to the weighted metric. The systematic use of the weighted metric and of the 
above orthogonal splitting as a fundamental scheme seems to be original. In particular, 
we show that the classical concepts of total kinetic energy, total kinetic momentum, total 
force, etc. can be regarded as a direct consequence of the above geometric scheme. 

Then, the formulation of a a constrained system of n particles can be easily obtained 
from the above scheme, by repeating the scheme of one free and constrained particle. 

Eventually, a particular care is devoted to the analysis of a system of n particles with 
a rigid constraint. 

First we study the geometry of the rigid configuration space, distinguishing the non 
degenerate and degenerate cases. Then, we formulate the kinematics and mechanics of 
a rigid system according to the above scheme. In particular, we show f Section 14. 2p that 
the classical formula of the velocity of a rigid system (well-known in mechanics) can 
be regarded as the parallelisation of a Lie group (well-known in differential geometry). 
This fact yields an interpretation of the inertia tensor as a representative of the weighted 
metric induced by the parallelisation. In this context, we exhibit a new explicit intrinsic 
expression of the angular velocity via the inertia tensor fCorollarv I4.24p . 

By combining the splitting of center of the mass and the splitting of the constrained 
configuration space, we obtain a splitting of the tangent and cotangent environmental con- 
figuration spaces into three components: the component of center of mass, the rotational 
component and a further orthogonal component to the configuration space (Theorem 
14.271 and Theorem I4.28|) . This splitting is refiected on all objects of the rigid system 
mechanics, providing a clear geometric interpretation of some classical constructions of 
mechanics and new results as well. For instance, the total momentum of forms arises from 
our geometric scheme via the projection on the rotational component fCorollarv I4.29p . 
Moreover, a special application of the above splitting is the explicit expression of the 
reaction force on every node fCorollarv l4.38|) . This formula seems to be new and possibly 
useful in engineering applications. 



Throughout the paper, we number those formulas which have a key role and/or em- 
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phasise a non standard feature of our approach. 

We assume all manifolds and maps to be C°° . If M and N are manifolds, then 
sheaf of local smooth maps M ^ iV is denoted by map(iW, iV) . 
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1 Preliminaries 

In this paper we use a few non-standard mathematical constructions. In order 
to make the paper self-contained, we start with some introductory notions. 

1.1 Scale spaces and units of measurement 

In order to describe in a rigorous mathematical way the units of measurements 
and the coupling scales, we introduce the notion of " scale space" fH] . 

We define a scale space U as "positive 1- dimensional semi- vector space" over 1R~^ . 
Roughly speaking, this has the same algebraic structure as , but no distinguished gen- 
erator over . We can naturally define the tensor product between scale spaces and or- 
dinary vector spaces. Moreover, we can naturally define the rational powers U^/^ of a scale 
space U . Rules analogous to those of real numbers hold for scale spaces; accordingly, we 
adopt analogous notation. In particular, we shall write 11° := U , U~"^ := U* , := . 

In our theory, these spaces will appear tensorialised with spacetime tensors. The scale 
spaces appearing in tensor products are not effected by differential operators, hence their 
elements can be treated as constants. 

A coupling scale is defined to be a scale factor needed for allowing the equality of two 
scaled objects and a unit of measurements is defined to be a basis of a scale space. 

We introduce the scale spaces T of time intervals, L of lengths and M of masses. 
We will consider time units of measurement mq G T , or their duals u*^ G T* . 

1.2 Generalised affine spaces 

Affine spaces are important for classical mechanics because they offer a geomet- 
rical model of the basic configuration spaces. Moreover, affine spaces constitute the 
appropriate framework for elementary differential calculus. In this paper, we need 
a more general definition of the standard notion. Namely, we introduce generalised 
affine spaces associated with (possibly non Abelian) groups. This generalisation is 
suitable for the description of the configuration space of rigid systems. 

A (left) generalised affine space is defined to be a triple {A, DA, I) , where A is a set, 
DA is a group and I : DA x A —>■ A is a. free and transitive left action. For the sake of 
simplicity, we often denote the generalised affine space {A, DA, I) just by A . 

For each a E A , the left translation lo '■ DA ^ A : g ^ go is invertible. In fact, for 
each a E A , there is a unique g E G , denoted hy g = ao^^ , such that a = go . 

A generalised affine map is defined to be a map f : A ^ A' between generalised affine 
spaces, such that, for a certain o E A ,we have /(a) = D/(ao~^)/(o) , for a\\a E A , where 
Df : DA —>■ DA' is a group morphism. We can easily prove that, if such a Df exists, then 
it is unique and independent of the choice of o . We say D/ to be the generalised derivative 
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of / . For example, if o G A , then the left translation : DA — > A is a generahsed affine 
map and its derivative is just the identity. 

Of course, if A is a generahsed affine space associated with the additive group of 
a vector space V , then A turns out to be an affine space according to the standard 
definition and the notions of generahsed affine map and generahsed derivative reduce to 
the standard ones. 

Now, let us consider a generalised affine space A associated with a Lie group G . Then, 
there is a unique smooth structure of A , such that the left translation / : G x A — A be 
smooth. 

Let Q := TeG be the Lie algebra of the group G , where e G G is the unit element. 
We recall the following well-known result (see, for instance, j2"5j). 

1.1 Lemma. The Lie group G is parallelisable through a natural isomorphism TG ~ 



Proof. Let us consider the map Tl : TG x TG TG and the trivial vector subbundles TeG C TG , 
over {e} C G , and G C TG , over G . 

Then, the restriction of Tl to the vector subbundle T^G x G yields, the linear fibred isomorphism 
TI\t^GxG -.TeGxG^TG over G . 

Thus, for each 5 G G , we have the linear isomorphism Tlglx^c : T^G TgG . QED 

We can easily generalise the above result to the affine space A . 

1.2 Lemma. The affine space A is parallelisable through a natural isomorphism 
~ A X . 

Proof. Let us consider the map Tl : TG x TA —^ TA and the trivial vector subbundles TeG C TG , 
over {e} C G , and AcTA, over A . 

Then, the restriction of Tl to the vector subbundle T^G x A yields, the linear fibred isomorphism 
TI\t,gxA ■■ TeG X A-^TA over A . 

Thus, for each a G A , we have the linear isomorphism Tlalr^c '■ TeG TaA . QED 

1.3 Lemma. Let M be a differential manifold equipped with a parallelisation TM ~ 
MxF , which induces a projection TTM — > TM . Then, we obtain the linear connection, 
given for each sections X, Y : M TM , by means of the composition 



G X g. 



M 



TM 



Y 



TM 



TX 



^ TTM .U 



Thus, the parallelisations of G and A induce linear connections of G and A . 
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2 Mechanics of one particle 

First, we review the one free and constrained particle mechanics as an introduc- 
tion to our formahsm and a pattern for next generaUsations. 

2.1 Free particle 

We start with a free particle moving in a 3-dimensional Euclidean affine space 
P . On the other hand, in many respects, the dimension 3 and the affine structure 
have no essential role; in fact, we essentially exploit the underlying weaker structure 
of Riemannian manifold of P . 

2.1.1 Configuration space 

We define the time to be a 1-dimensionaI affine space T associated with the vector 
space T := T (g) iR. We shall always refer to an affine chart (x") induced by an origin 
to & T and a time unit of measurement m° G T . 

We define the pattern configuration space to be a 3-dimensional affine space P asso- 
ciated with an oriented vector space S . We shall refer to a (local) chart (x*) on P . Latin 
indices i ,j ,h ,k will run from 1 to 3 . 

We shall also be involved with the tangent space TP = P x S and the cotangent 
space T*P = Px S* . We shall refer to the local charts {x\ x') of TP and {x\ Xi) of T*P 
and to the corresponding local bases of vector fields di and forms . We also denote by 
{x\ x\ x\ x^) the induced chart of TTP , with {di, di) and {d\ d^) the corresponding bases 
of vector fields and 1-forms; we have the chart of VTP . 

The parallelisation of P induced by the affine structure yields a flat linear connection 
V (see Lemma [Ol) 

We equip S with a scaled Euclidean metric g ® {S* <S) S*) , called pattern metric, 
which can be regarded as a scaled Riemannian metric of P 

g-.P^V® {T*P ® T*P) . 

We denote by g the corresponding contravariant metric. We have the coordinate 
expressions g = gijd' ® d^ and g = g'^^ di ® dj , with gij G map(P, ® M) and 
g^^ G map(P, (g) M) . The associated flat isomorphism and its inverse, the sharp 
isomorphism, are denoted by g^ : TP — > T*P and g'^ : T*P TP . The 

metric g and an orientation of S yield the scaled volume form t] G L'^ (g) A^S* and its 
inverse fj G L*^ (g) . 

The Riemannian connection associated with g coincides with V . We denote the ver- 
tical projection associated with V by z/ : TTP TP and the the Christoffel symbols by 

■pi 

hk ■ 



2.1 Free particle 
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2.1.2 Kinematics 

We define the phase space as the 1st jet space of maps T ^ P 

JiP := T X (T-^ ® TP) = TxPx (T"^ ® S) . 

The induced chart of JiP 

A motion is defined to be a map s : T P . 

The 1st differential, the 2nd differential, the velocity and the acceleration of a motion 
s are defined to be, respectively, the maps 

ds:T^T-^0TP, d^s:T^ T'^ ® r(T-^ ® TP) , jis J^P, 

Vds : T ^ TP . 

By definition, we have jis(t) = (t, ds(t)) and Vds = v o d^s . 
Moreover, by taking into account the splittings 

TP ~ P X (T~^ ® S) , 
T-^ ® T{T-^ ® TP) ~ (P X (T-^ ® S)) X ((T-^ ® S) x (T'^ ® 5)) , 

we can write 

ds = {s, Ds) , dh = {s, Ds, Ds, D^s) , Vds = {s, D^s) , 

where Ds : T — > T"-*^ 5 is the standard derivative of s . 
We have the coordinate expressions 

{x\xq) ods = {s\ dos') , 

{x\xi,xi,xio) ° = 9os\ dos\ d^s') , 
{x^,x\x'o)ojis^ s\ dos') , 
{x\ 4o) o Vds = {s\ d^s' + (4, o s) dos^ dos'^) . 



With reference to a mass m e M , we define the kinetic energy and the kinetic mo- 
mentum, respectively, to be the maps 

: T"^ ® TP ^ ®M) ^ M : v ^ \ m g{v, v) 

V := DIC : J-^ (g) TP ^ (T"^ (g O M) O T*P : v ^ m g\v) . 

We have the coordinate expressions K, = \ m gij x^ Xq and V — m gij Xq d^ . 



12 



2 Mechanics of one particle 



2.1.3 Dynamics 

In our context, the force acting on a particle is given a priori on the the phase space. 
Moreover, it is convenient to introduce the force as a co-vector. 
Thus, a force is defined to be a map 

F : JiP ^ (T-^ (8) (8) M) (8) T*P . 

The force F is said to be conservative if it factorises through JiP ^ P and can be 
derived from a potential U : P ^ M iR by the equahty F = dU .li the force 
is conservative, then we define the associated Lagrangian to be the map C :— K, + U . 

We say that a motion s fulfills the Newton's law of motion if 

(2.1) mg^{Vds) = Fojis. 

It is remarkable that we can link the formulation of dynamics in terms of the con- 
nection V with the Lagrangian approach, directly without any reference to variational or 
Lagrangian calculus. In fact, the following Lagrange's formula holds 

m g^ (Vds) = {p{diK. o ds) - dilC o ds) {d' o s) . 

By the way, the above formula provides quickly the Christoffel's symbols of V . 
Hence, the coordinate expression of the Newton's law is 

m gij (a^s^' + {Ti^ o s) dos'' dos'') = D{diK. o ds) - diJC ods^FiO j^s . 

In the particular case when the force is conservative, the Newton's law of motion is 
expressed by the Lagrange equations 

D{diC o ds) — diC o ds — . 
2.2 Constrained particle 

We assume an embedded submanifold of the pattern Euclidean affine space as 
configuration space of a constrained particle. 

The mechanics of a constrained particle has two features: an intrinsic and an 
extrinsic one. According to the intrinsic viewpoint, the particle behaves as a 'free' 
particle moving in an Z-dimensional Riemannian manifold; hence, according to the 
intrinsic viewpoint, we can repeat the scheme of the previous section. On the other 
hand, the environment space adds an exterior geometric structure: the 2nd funda- 
mental form, which measures the deviation of the sTibmanifoId from being an affine 
subspace of the environmental space. Then, according to the extrinsic viewpoint, we 
interpret the reaction force in terms of the 2nd fundamental form of the constrained 
space. 



2.2 Constrained particle 
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2.2.1 Configuration space 

We define the configuration space for a constrained particle to be an embedded sub- 
manifold of dimension 1 < / < 3 

, ■ p ^ P 

''con ■ con ' ■ 

Thus, by definition of embedded submanifold, for each p G -Peon , there exists a chart 
(a;*) of P in a neighbourhood oip , such that Peon is locally characterised by the constraint 
l^z+i _ ,j,3 _ Qj. ^jigji^ (^yi^ ■■■,?/') — {x^\Pcon: • • • ) ^'|Pcon) tums out to be a local 

chart of Peon ■ The functions y^, . . . ,y'' are said to be local Lagrangian coordinates and 
the functions x'^'^^, . . . to be local constraints. Prom now on, we shall refer to such 
adapted charts. 

For practical reasons, we shall adopt the following convention: 

- indices i, j, h, k will run from 1 to 3 ; 

- indices a,b,c,d will run from 1 to Z ; 

- indices r, s, t will run from / to 3 . 

We have TP\p_ = Peon x S and T*P|p^„„ = Peon X 5* . 
We have the natural dual linear injective and projective maps 

Tieon : TPeon ^Pjp^ C TP and TT := TXon : r*P|p_ ^ T*Peon , 

with coordinate expressions 

rieon( E X"da)= ^ ^" ^nd 7l{ ^ '^i^')^ E ^ad\ 

l<a<l l<a<l l<i<n l<a<l 

The complementary linear projective and injective maps can be obtained by means of 
the metric. 

We consider the orthogonal subspaces 

T^Peon {X e rP|p_ I g{X, TPeon) = 0} C TPjp^ , 
TxPcon := {a e r*P|p_ I a(rPeon) = 0} C T*P|p_ . 

The vector fields da are tangent to Peon , while the vector fields dr are transversal. If 
dr e T^P con , then the adapted chart (x*) is said to be orthogonal to the submanifold. 

The subspace T^Pcon consists of the forms of the type cu — J2i+i<r<n^rd^ , i.e. of 
forms whose "tangent" components vanish. 

The restriction 

^con := i*con9 ■ ^con ^ ® (r*P,on « T*P,on) 

of the pattern metric g to Peon is a scaled Riemannian metric, which will be called the 
intrinsic metric. Its coordinate expression is gcon = J2i<a,b<i{9con)abd"' ® d'^ , where we 
have set (fi'con)o6 •= 9ab\Pcon ■ The contravariant form of g^con will be denoted by ^con • We 
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stress that, in general, the / x / "tangent" submatrix of (g^^) is different from the inverse 
of the matrix {gab) ; they are equal if and only if the adapted chart (x') is orthogonal. 
The intrinsic metric gcon yields the Riemannian connection Vcon ■ 
With reference to a mass m e M , we define the intrinsic kinetic energy and the 
intrinsic kinetic momentum 



Peon := i*V : T-^ ® TP, 



R : V ^mgcon{v,v) , 
•T*P : V mg\on{v) , 



with coordinate expressions /Ccon = | g^on ah v'' and Vcon = m g^on ab V"" d'' . 

Let us analyse the orthogonal splittings of the tangent and cotangent spaces induced 
by the metric. 

The metric yields the injective map 



J : T*Pc 

through the commutative diagram 



T* P\ 



T* P\ 



3' con 



L- 



L- 



T*P\ 



We have the mutually dual orthogonal splittings 



TP|p_ = TP 
T*P\p = T*P, 

■*■ con ' 



con © T Peon 7 

-f*con 



TlPr 



with projections 



TP 

con ) 

rp* p 



con ) 



7r± : T*P\p^ 



T^Pc 

TlPr. 



As the projection tt has a very simple expression, it is convenient to compute the other 
projections tt", tt"*-, tt± via the following commutative diagrams 



TP I 



L2 (g) T*P| 



TP 

J- r. 



9' con 



L2 ® T*Pp 



2.2 Constrained particle 
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7lt 



9^' 



Then, for each X e TP|p^^^ and a; e T'*P|p^^^ , we obtain the equahties 

l<a,b,c<l 

^"W= E E X^^.6<7ct5a= E E {X'^ + X^gr.g'^^Jda, 

l<i<n l<a,b<l l<a,b<l l+l<r<n 

^HX)= E E X^idr -9rb9'.Lda), 

l<a,b<l l+n<r<l 

7r± (t^) = E E (^^ ~ ^» ^con 9br ) d'' . 

l<a,b<l l+l<r<n 

Of course, the above formulas simphfy considerably if the adapted chart is orthogonal, 
i.e. if ^rblPcon = 0. 

2.2.2 Kinematics 

We define the intrinsic phase space as the 1st jet space of maps T — > Peon 

JlPeon := T X ® TPeon) • 

The induced chart of JiPcon is y", yo) ■ 

A constrained motion is defined to be a map Scon : T Peon C P . Clearly, a 

constrained motion can be naturally regarded as a motion of the pattern space, via the 
inclusion icon • Indeed, a motion s : T — > P is constrained if and only if s** = . 

The 1st differential, the 2nd differential and the velocity of a constrained motion 
Scon , computed in the environment space, turn out to be valued in the corresponding 
constrained subspaces 



dSeon : T ^ TPeon , C^'Scon : T ^ T'^ ® T{T-' ® TPeon) , 

Jl^con ■ T > JlPcon • 

We have the coordinate expressions 

(l/^2/o)°a?Scon = (s^aoS»), 

(2/^^o^i/S,i/oo) ° c^^^eon = 5os^ 9o.^ ^o^^'^) 

(a;^|/^|/o")oJlSeon=(x^ 5», ^oS'^). 

Hence, as far as the above objects are considered, the intrinsic and the extrinsic 
approaches coincide, up to the natural inclusion of the constrained spaces into the corre- 
sponding environmental spaces. 
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Conversely, the intrinsic and extrinsic approaches of the acceleration of the constrained 
motion Scon do not coincide. The intrinsic viewpoint is suitable for the intrinsic expression 
of the law of motion and the extrinsic viewpoint provides the constraint reaction force. 

We define the intrinsic acceleration of a constrained motion Scon as the map 

with coordinate expression 

VconC^Scon = {^0^'' + (^con tc) ° ^con ^O^^ ^O^'') (<9a O Scon) ■ 

l<a,b,c<l 

Analogously to the free case, the intrinsic co-acceleration is given by the Lagrange's 
formula 

mgconi^conds) = {p{daK,con ° dScon) - <9a/Ccon O dScon) {d" O Scon) • 

On the other hand, by regarding the constrained motion Scon as a motion of the 
environmental space, we define the extrinsic acceleration as the map 

Vdscon : T ^ O TP , 

with coordinate expression 

l<j<3 

VrfScon = (^0*' + O^bc) ° Scon ^O*^ ^QS") (di O Scon) • 

l<a,b,c<l 

Then, according to the Gauss' Theorem [7j, we have the splitting 

VfiScon = Vr"(V(i,Scon) + vr-^(V(i,Scon) , 

with 

7rll(VdScon) = VconfiScon and TT-^ O VrfScon = O dScon , 

where 

N : TP|p_ ^ T^Pcon 
is a quadratic map, called 2nd fundamental form, whose coordinate expression is 

N= Y E Tlj'r{dr-9rb9'.M. 

l<a,b,c<l «+l<r<3 

Indeed, the map N measures how the submanifold Peon deviates, at 1st order, from 
being an affine subspace of P . The quickest way to compute the 2nd fundamental form 
is the following: compute the covariant expressions of the extrinsic and intrinsic acce- 
lerations via the Lagrange's formulas; then pass to the contravariant expressions and take 
the difference. 



2.2 Constrained particle 
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2.2.3 Dynamics 

Let us consider a force F : JiP — >• (T~^(8)L^(8)M) (g)T'*P in the environment space P . 
As we are dealing with constrained mechanics, we are involved only with its restriction 

F ^1 j,p_ : JiPeon ^ ® ® M) ® r*P|p_ . 

According to the splitting of T*P\p^^^ , we can write 

F — F + F 

^ — ^ con \ ^ con 1. i 

where 



Feon = i*F : JiPeon ^ (T"^ M) T*P^^^ . 

We call Fcon the intrinsic force. 

Let us assume that constraint confines the motion on the configuration space Peon , 
via Newton's law of motion, by means of a suitable additional 'reaction force' defined on 
the constrained space 

R : JiPeon ^ (T~2 ® ® M) ® T*P\p^^^ . 
According to the sphtting of T*P\p^^^ , we can write 

where 

i?con = i*R : JlPcon ^ {T'^ 0) M) T*Peon ■ 

We call i?con the intrinsic reaction force. 

A constrained motion Scon : T — > Peon is said to fulfill the constrained Newton's law 
of motion if the following equation holds 

mg^ O (VdScon) = {F + R)o jiScon • 

2.1 Theorem. A constrained motion s : T — > Peon fulfills the constrained Newton's 
law of motion if and only if 

"^5'con ° (VconC^Scon) = (-^con + -^con) ° Jl^con 

mg^ o{N O dScon) = {Fcon ± + ^con ±) ° jl^con • 

Actually, for each choice of initial data in JiPeon , the 1st equation has (locally) a 
unique solution and the 2nd equation is fulfilled if and only if 

(2.2) Rcon±^mg^ oN - F^on±-a 
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2 Mechanics of one particle 



According to the above result, we make the "minimal assumption" {virtual works 
principle, i.e., smooth constraint) 

Rcon ■ 

Then, the explicit coordinate expression of R is 

l+l<r<3 

(2.3) = ^ (m {Tarb - grc 9con Tcon adb) 2/" - Fr + Q^l^ Fa Qbr) O icon (c^" O . 

l<a,b,c<l 

In classical literature (see, for instance, [13]) the computation of reaction force is pre- 
sented implicitly as the solution of a linear system associated with Lagrange multipliers. 
Instead, the above formula (which involves an adapted chart) provides an explicit ex- 
pression of the reaction force in terms of the Christoffel symbols, or of the metric. In the 
particular case when the adapted chart is orthogonal, this formula becomes very easy. 

Thus, the dynamics of a constrained particle can be interpreted as the dynamics of 
a 'free' particle moving on the Riemannian configuration space Peon • Then, all main 
notions and results holding in the free case by using the Riemannian structure of P can 
be easily rephrased in the constrained case. This is a remarkable conceptual and practical 
advantage of the present approach. 
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3 Mechanics of a system of n particles 

In this section we generalise the previous concepts and results to systems of many 
particles. Our guideline will be the interpretation of the multi-particle system as a 
one-particle moving in a higher dimensional space. In this way, all we have learned 
for one-particle can be applied directly to multi-particle systems. On the other hand, 
we have additional concepts, e.g. the center of mass splitting, which follow from the 
projections on the factor spaces of the different particles. 

3.1 Free particles 

We start with a finite system of free particles. 

The key feature of our analysis is the geometry of the configuration space, which 
is naturally equipped with two metrics and two natural splittings. 

We shall use systematically the prefix "multi" to indicate objects of the n-system 
analogous to objects of one-particle. 

We assume n > 1 , and consider n masses mi, ... , m„ G M . 

For each i — 1, . . . ,n , with reference to the i-ih particle, it is convenient to consider 
a copy of the following pattern objects: = P and Si = S . 

3.1.1 Geometry of the multi— configuration space 

The multi- configuration space is defined to be the product space 

-Pmul := Pi X ■ ■ • X P„. 

Clearly, Pmui is an affine space associated with the vector space Smui '■— -Si x • • • x 5„ . 

A product chart (x^^i) of Pmui induced by charts {x\) , . . . , (x^) of the single com- 
ponents is said to be without interference of the particles. Conversely, a chart {xl^^^i) of 
Pmul which cannot be written as a product as above is said to be with interference of the 
particles. 

We define the total mass as mo X^^mj e M and the i-th weight as //j = mj/mo . 
Clearly, we have J2i fJ'i = ^ ■ 

A typical notation for the elements of Pmui , -S'mui and S^^^ will be 

Pmul = (Pl, ■ ■ ■ ,Pn) e Pmul = Pi X 
Vmul = {Vl, ■■■,Vn) £ 'S'mul ^ Si X 

ttmui = (oil, . . . , «„) e 5^ui = 5t X 



... X P„ , 

... X Sji , 
... X 5! . 
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3 Mechanics of a system of n particles 



We define the multi- geometrical metric and the multi-weighted metric as 

(3.4) g^nnl ■ Pmnl ^ ® (T*Pmul ® T*Pmul) : {Umnh ^mul) ^ ^jdi^U Vi) ; 

i 

(3.5) Ginul : -Pmul ^ ® (T*Pmul ® 7'*Pmul) : (^imul, ^'mul) ^ J2 l^i di^i, Vi) . 

i 

The contravariant tensors of (jfjnui and Gmui are denoted, respectivefy, by ^mui and Gmui ■ 
If = 1 , then Gmui = fi'mui = 9 and m = rriQ . 

If n > 2 , then the two metrics Gmui and gm_vi are distinct. Moreover, if mi = ■ ■ ■ = rrin , 
then = = = i and Gmui = ^ fi'mui • 

For a system of n particles, we will rephrase the dynamics of a system of one particle, 
by replacing the pattern metric g with the weighted metric Gmui and the mass m with the 
total mass mo . This procedure yields the correct Newton's law of motion, in full analogy 
with the one particle case. 

According to this scheme, we define the multi-kinetic energy and the multi-kinetic 
momentum by 

^mui : ® TP^ui ^ (T-^ ®V(^m)®IR : ^ i mo G^^i{v, v) , 
Vrani : ® TP^^i {T'^ ®l} ®M) ® T*P^^i -.v^nio G^^^iiv) , 

and recover the standard formulas /Cmui(^') = J2i ^m-i g{vi,Vi) and Vinu\{v) = (mig{vi)^ . 

The linear connection Vmui induced on the multi-configuration space Pmui by the affine 
structure (see Lemma II. 3p coincides with the Riemannian connection induced by both 
metrics fi'mui and Gmui • This is true although the two metrics need not to be proportional. 



3.1.2 Splittings of the configuration space 

The multi-configuration space has two distinguished splittings. 

Multi— splitting. We have the obvious affine multi- splitting Pmui = Pi x . . . x P„ . 
The corresponding affine projections vr, : Pmui Pi and the further induced projections 
can be used to extract information on the single particles from the kinematical and 
dynamical multi-objects of the multi-system. 

The subspaces Si, . . . , Sn C S^ui are mutually orthogonal with respect to both met- 
rics fifmui and Gmui • 

Diagonal splitting. Moreover, Pmui has the following further diagonal splitting, 
which has no analogous in the one-particle scheme. 

We define the diagonal affine subspace, the diagonal vector subspace and the relative 
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vector subspace, respectively, as 

-Pdia '■— {Pmul ^ -Pmul | Pi = ' ' • = Pn} C Pmul , 
'S'dia ■= {'f^mul £ 'S'mul \ Vi — ■ ■ ■ — Vn} C Smul , 

i 

In the following, the subscripts "dia" and "rei" will denote the objects associated with 
the above spaces. 

3.1 Theorem. We have the affine splitting of the multi-configuration space and the 
linear splittings of the associated vector and covector multi-spaces 

(3.6) Pmul = Pdia © S,el ■ Pmul = (Po, • • • , Po) + (Pl " Po, • • • , Pn " Po) , 

(3.7) S'mul = S^dla © -S^rel : 't'mul = {vq, . . . , Vq) + {vi - Vq, . . . , Vn - Vq) , 

(3.8) = Sli^ © : a^nnl = (/^l ^O, ■ ■ ■ , l^n ^o) + ("i - /^i ao, /^n tto) , 

where, for each Pmui £ Pmui , Po & P is the unique point such that 

i 

and where, for each Vmui £ -S^mui o.n'd amui £ -S^mui > 

i i 

The above splittings are orthogonal with respect to the weighted metric Gmui ■ 
Proof. 1) Let us prove the splitting S'mui = 'S'dia © Srei . 

For each vq G S , we have J2i A*«'^o = if and only if i>o = . Hence, Saia H Srei = . 

Moreover, for each Vmui S S^mui , we have {vi, . . . , v„) = (wq, ■ • • , vq) + (vi —vq, . . . , w„ -Vq) , with any 
Vo G S . Clearly, {vo, . . . ,vo) € S'dia and {vi—vo — vq) G Siei , if and only if vo = J2i l^i'^i ■ Hence, 

-Sdia + Sre\ = Smui and the expression of the splitting is expressed by the 2nd formula of the statement. 

2) Let us prove the splitting Pmui = -Pdia © Srci ■ 

Clearly, P^ia C Pmui is an afhne subspace associated with the vector subspace S'dia C Smui • Hence, 
the equality Smui = S^ia,® Sre\ impUes Pmui = -Pdia ®S,ei- 

Moreover, for each p„,^i G P^ui , we have (pi, . . . ,pn) = {po, ■ ■ ■ ,Po) + {Pi ~Po, ■ ■ ■ ,Pn-Po) , with any 
Po ^ P - Clearly, (po, • • • ,Po) G -Pdia and {pi - Po, ■ ■ ■ ,Pn - Po) G "^rei if and only if J2i l^i {Pi - Po) = . 
Hence, the expression of the splitting is expressed by the 1st formula of the statement. 

3) The splitting ^mui = S'dia © S^d impHcs the splitting S*^^^^ = ^dia © S*^^ . 

Moreover, for each Umui G S'^^i , we have {ai, . . . , a„) = (/3i, . . . , /?„) + (ai - /3i, . . . , a„ - /?„) , with 
any A e S* . 

On the other hand, for each v^^ui € -Smui , we obtain 

(/3i,..., /?„)(«!,..., t;„) = {ai,...,an){vo,...,vo) 

if and only if 

i j i i 3 
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i.e., if and only if 

3 

Moreover, in virtue of the equalities 

i j i j i i 

we obtain 

(ai - /i^i (X"i)' ■ • ■ ,an - C^aj)){vi,. ..,«„) = (ai, . . . ,a„)(ui - uq, . . . , u„ - uo) • 
j j 

Hence, the expression of the splitting is expressed by the 3rd formula of the statement. QED 
We denote the projection associated with the above sphttings by 



TTdia : 


■ .^'mul 


-f dia ) 


TTrel : 


■ .^'mul 


'S'rel 


TTdia : 


■ 'S'mul ~ 


"S^dia ) 


TTrel '■ 


■ 'S'mul 


'S'rel 


=.* 

^dia • 


■ '-'mul ~ 


"^dia > 


=.* 
rel ' 


■ '-'mul 





Clearly, the above splittings depend on the choice of the multi-mass and are not 
orthogonal with respect to the geometrical metric (unless all masses are equal) . 

We stress that, while we have the natural inclusion Pdia ^ -F*mui , we do not have a 
natural inclusion S'rei ^ Pmui ■ 

We have a natural splitting of the weighted multi-metric of the type Gmui = GdiaBG^rei • 
Moreover, the affine structures of Pdia and Sj-^i yield the flat connections Vdia and Vrei , 
which turn out to be the Riemannian connections induced by Gdia and Grei > respectively. 

Center of mass splitting. We can describe the diagonal splitting in another way, 
via the center of mass. 

According to the above Theorem, we define the center of mass of pmui £ -Pmui to be 
the unique point Pq ^ P , such that J2i f^i{Pi ~ Po) = . By considering any o e P , wc 
can write Po = o + fJ'i [Pi ~ o) . With reference to the center of mass, it is convenient to 
consider a copy of the following pattern objects: Pcen = P , -S'cen = S and S^^^ = S* . 

Thus, we have the center of mass affine projection 

TTcen : -Pmul ^ -Pcen : Pmul ^ Po O + ^ fXi {pi - o) , for any O G P . 

i 

The linear projections associated with the afline projection TTcen ■ -Pmui ~^ Pcen turn 
out to be, respectively, the weighted sum and the sum 

i 

Seen : S*^^^ Sl^^ : Qimui ^ ao ■.= Y,oti- 
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The 2nd projection is just the map which associates with every multi-form its total 
value. Indeed, this map plays an important role in mechanics of multi-systems. 
Clearly, we have the natural affine isomorphism and linear isomorphisms 

Pcen fdia : Po ^ (Po, ■ ■ ■ , Po) , 
Seen 'S'dia : ^^0 ^ {Vq, ■ ■ ■ , Vq) , 

3.2 Corollary. We have the center of mass splittings 

Pmul - Pcen X Src\ ■ Pmul - [pO , (Pl - PO, • • • , Pn " Po)) , 
S'mul ^ Seen X S^el ■ t^mul ^ (^^0 , {Vi - Vq, . . . , Vn - Vf))) , 
-^mul - -^cen ^ ^*e\ ■ "mul - ("O , ("l " yUl "O, • • • , " /^n "o)) ■ □ 

According to our scheme, we define the center of mass kinetic energy and the center 
of mass kinetic momentum as 

JCeen ■ ® TP^en (T^^ (g) V (g) M) (g) R : Veen ^ ^ mQ QeeniVeen, Veen) , 

Veen ■ ® TPeen ^ (T^^ ® ® M) ® T*Peen ■ Veen ^ QeeniVeen) ■ 

Moreover, we obtain the equalities 

^mul(^^mul) = ^cen{Vo) + /Cmul ("^^rel) , 

JCeen ° T^een = f^mul ° ^^TTdia : ?^mul ^ \ rUQ QeenivQ, Vq) = \mQ^lli g{Vo, Vq) , 

i 

Veen ° THeen = Seen ° ^mul ° ^^TTdia : ^mul ^ "^0 fi'cenl^o) = "^0 XI 5'^(^o) ■ 

i 

3.1.3 Kinematics 

For the kinematics of our system of n particles, we follow the viewpoint of a multi- 
particle moving in a multi-space (in analogy with a one-particle moving in the standard 
space). Moreover, we take into account the center of mass splitting. 

We define the multi-phase space as JiPmui '■= T x (T^^ ® TPmui) • 
A multi-motion is defined to be a map Smui : T — * Pmui • 

Of course, the multi-motion can be regarded as the family of motions of the system: 
Smui = (si, Sn) ■ This holds also for the derived quantities, like the multi-velocity 
dsmui : T (g) TPjnui and the multi- acceleration VmuiC^Sinui : T T^^ eg) TPmui • 

We can relate the multi-motion to the splitting of center of mass (Corollary 13. 2|) by 
the equalities 

"Smul "^dia ~l~ '^rcl — ("^cen y "^rel) ; ds^nul '^'5^;^ ~l~ ds^el — i^dSeen ; '^■^rel)- 

^mul'^'^mul ^dia '^■^dia ^rel dSj-gi — (Veen dSeen ; ^rel '^■^rel) • 
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3 Mechanics of a system of n particles 



3.1.4 Dynamics 

In analogy with the case of one-particle, we define a multi-force to be a map 

i^mul : JlPmnl ^ (T"' (g) ® M) ® T^P^^i . 

Again, the multi-force can be regarded as the family of forces acting on each particle 
-fmui = (-^1, • • • , Fn) . In general, each of the components is defined on the whole phase 
space. In the particular case when each component Fj of the multi-force depends only the 
i-th phase space, the multi-force is said to be without interaction. 

We say that a multi-force Fmui fulfills the Newton's 3rd principle if, for each 1 < i < n , 

Fi= ^ Fij , Fij{pi,pj) = Xij{\\pj - Pi\\g) g\pj - Pi) , Xij = Xji , , 

l<i^j<n 

where Fij : Pi x Pj {T'^ O O M) O T*Pi and Xij : R ^ IR , for each 

1 < i ^ j < n . 

The total force of the system is defined to be the component of the multi-force with 
respect to the center of mass 

i 

The multi-force is said to be conservative if it can be derived from a multi-potential 
Wmui '■ -Prnui — (g) (g) M) (8) 1? as Fmui = dUinui ■ this case, we define the multi- 
Lagrangian to be the map 

We say that a multi-motion Smui fulfills the Newton's law of motion if 

(3.9) mo G^^i O (Vmul ds^^i) = Fmul O JlSmul ■ 

We can split the Newton's law with respect to the multi-splitting and to the splitting 
of the center of mass. 

In the former case, we simply obtain the system of coupled equations 

rrii g^ o dsi = Fj o jis . 

In the latter case, we have the following Theorem. 

3.3 Theorem. The Newton's equation is equivalent to the system 

(3.10) mo fi-cen ° (Veen ds^en) = -^cen O Jl^mul , rUQ G^^^ O (Vrel ds^el) = F^^i O jiS^ul ■ □ 

If Fcen factors through JiPcen , then the 1st equation can be integrated independently 
of the 2nd one and can be interpreted as the equation of motion of the center of mass. 
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As for one-particle, the coordinate expression of the Newton's law is, with reference 
to any chart (x^^j) of Pmui , 



which, if Fjnui is conservative, is equivalent to the system of Lagrange's equations 



3.2 Constrained particles 

According to our programme, the analysis of the geometry of the constrained 
space for a system of n particles can be carried out by analogy with the case of 
one-particle. 

We assume the multi- configuration space of a constrained system of n particles to be 
an embedded submanifold Pmui con C Pmui • 

In general, it is not possible to write Pmul con = -Peon l X . . . X Peon n , with Peon i C Pj , 

for each i = 1, . . . ,n . In the particular case when this holds, we say that the constraint 
is without interference between particles. 

Moreover, in general, it is not possible to write Pmui con = Pcen con x -S'rei con , with 
Peon con C Pcen and S^ci con C S^ci ■ In the particular case when this holds, we say that 
the constraint is without interference between center of mass and relative positions. In 
this case, the intrinsic metric Gcon splits into the sum of the metrics Gcen con and Grei eon 
(according to Corollarv 13. 2|) . with interesting consequences in dynamics. 

We leave to the reader the task to formulate the kinematics and dynamics of a con- 
strained system of n particles according to our scheme. 
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4 Rigid systems 

Now, we specialise the theory of constrained systems of n particles to the case 
of a rigid constraint. We devote emphasis to the geometric structme of the rigid 
configuration space. In particular, we show that this is the true source of the classical 
formulas of the velocity of rigid systems. 

Throughout the section, we suppose that n > 2 . 

4.1 Geometry of the configuration space 

Let us define the scaled functions, for 1 < i, j < n , 

Tij ■■ Pmul ^ L (g) iR : Pmul = (Pl, • • • ,Pn) \\Pi " PjWg ■ 

A rigid configuration space is defined to be a subset of the type 

Vig • -frig • {Pmul ^ mul | ^ij{Pm.u\) t ^ — ^ ^ j — ^} C -Pmul ) 

where e L fulfill 

li3 = lji, l<i:j<n, i^j, 

kk < kj + ljk, I <i,j,k <n, iy^j,jy^k,kj^i. 

Note that we have excluded the case in which the positions of different particles 
coincide. 

Prom now on, let us consider a given rigid configuration space Prig . 
We define the rotational space to be the subset 

iiot '■ Srot — {Vrel £ 'S'rel | l^i " ^ kj , 1 < 1 < j < u} C. S^el ■ 

A typical notation for the elements of S^ot will be 

^rot = (ri, ...,rn) e 5rot C S'rel ■ 

Due to the equality \\pi — Pj\\ = ||7rrei(Pmui)i — '?'"rei(Pmui)j|| = hj , the rigid constraint 
does not involve the center of mass but only relative positions. Then, the restrictions of 
the projections tTcch : -Pmui -Pccn and TTrei : -Pmui 'S'rei to the subset Prig C Pmui turn 
out to be, respectively, projections 

^cen ■ Prig ^ Pcen and TTrot ■ -^rig ^ ^lot • 

Thus, we obtain the bijection 

(^cen ) ^rot) ■ -^rig ^ Pcen ^ ^rot • 

Next, we classify the rigid constraints as follows. 



4.1 Geometry of the configuration space 
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For each t>rei G S^c\ , we define the vector subspace 



(■^rei) •= span{t>j — Vj I I < i, j < n} C S . 



4.1 Lemma. For each v^ei G S^d , we have (frei) = spanl^j — Vh \ h = 1, . . . ,n} , ior 
any chosen 1 < i < n . More precisely, Vi = J^hj^i f^h {vi — Vh) ■ 



If Trot G S'rot , then we call {r^ot) the characteristic space of (rrot) • 
Let us prove that the dimension of the characteristic spaces does not depend on 
elements in Srot , but only on Srot • We need the following technical Lemma (see also 0). 

4.2 Lemma. Let Trot , '"rot ^ "^rot ■ Then, there is an isometry (p : S ^ S such that 
0(rj) = r[ , for each i = 1, . . . , n . 

Proof. Let us fix i G {1, . . . , n} , and set dim (riot) — I ■ We can choose a basis of (wrot) of the type 
{hh) ■— {ri — \ h = 1, . . . ,1} . Consider the subset (6'^) := {r[ — r[^^ \ h = 1, . . . , Z} . We can define a 
linear map : {r^ot) (Kot) by 4'{bi) = K i ^'^^ ea.ch h = 1, . . . ,1 . We have g{(f>{bi),(p{bj)) = g{bi, bj) , 
in virtue of the equalities \\bh\\ = hh — \Wh\\ and \\bh — bk\\ — Ihk = — b'^W ■ The map is an 
isometry between (rj-ot) and {r[^^.) , hence it can be extended to a linear isometry of S . Eventually, we 
have g((/){rj), 4>(bh)) = g{rj, bh) , and the rigid constraint yields g{rj, b^) = 4>{i>h)) ■ Hence, we obtain 
4'iTi) = '^'j J for each j = 1, . . . , n . QED 

Therefore, if r^t , ^^t ^ -^rot , then dim {r^ot) = dim (r^^J . 

We define the characteristic of Prig to be the integer number Cp^^^ := dim (rrot) , where 
Trot £ S'rot • Obviously, we have 1 < Cp^.^ < 3 , and we can classify the rigid configuration 
space in terms of Cp^.^ . 

We say Prig to be 

- strongly non degenerate if Cp^.^^ = 3 , 

- weakly non degenerate if Cp^.^ = 2 , 

- degenerate if Cp . = 1 . 

^ J- rig 

Of course, if n = 2 , then Prig is degenerate; if n = 3 , then Prig can be degenerate or 
weakly non degenerate. 

Thus, by considering all particles as assuming positions in the same space P , the above 
cases correspond respectively to the case when the minimal affine subspace containing all 
particles is a line, or a plane, or the whole P. As a consequence of the above result, the 
case occurring for a given rigid system does not change during the motion. 

Let us denote by 0{S,g) the group of orthogonal transformations of S with respect 
to g . We want to study the topological subspace Srot C Srei through the natural action 
of the Lie group 0{S,g) on Srot • More precisely, we can easily prove that the map 



Proof. We have -fiiVi = Y^h^i f^^Vh and Vi ~ fiiVi 
tion. QED 



J2h=ii i^hVi , which give the result by subtrac- 
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is well-defined and yields a continuous action of 0{S, g) on S^ot ■ Such an action of 0{S, g) 
on Srot is transitive because of Lemma 14.21 Let us denote the isotropy group at rmt by 
H{T,,,)(lO{S,g). 

4.3 Lemma. The following facts hold: 

- in the strongly non degenerate case the isotropy subgroup H[r.^at\ is the trivial sub- 
group {1} ; 

- in the weakly non degenerate case the isotropy subgroup H[r ^ot] is the discrete sub- 
group of refiections with respect to (rrot) ; 

- in the degenerate case the isotropy subgroup H[r-[ot] is the 1 dimensional subgroup 
of rotations whose axis is (rrot) ; we stress that this subgroup is not normal. □ 

4.4 Proposition. The following facts hold: 

- Srot is strongly non degenerate if and only if the action of 0{S,g) on ^rot is free; 

- Sj-ot is weakly non degenerate if and only if the action of 0{S, g) on S^ot is not free, 
but the action of the subgroup SO{S,g) C 0{S,g) on S^ot is free; 

- Srot is degenerate if and only if the action of SO{S,g) on S^t is not free. □ 

Of course, if n = 2 , then S^ot is degenerate; if n = 3 , then S'rot can be degenerate or 
weakly non degenerate. 

Furthermore, we can prove the following result (recall the definition of affine space in 
Section ir^ . 

4.5 Corollary. The following facts hold: 

- if Srot is strongly non degenerate, then Srot is an affine space associated with the 
group 0{S,g) ; 

- if Srot is weakly non degenerate, then S^ot is an affine space associated with the 
group SO{S,g); 

- if ^rot is degenerate, then S^ot is a homogeneous space (i.e. the quotient of a Lie 
group with respect to a closed subgroup) with two possible distinguished diffeomorphisms 
(depending on a chosen orientation on the straight line of the rigid system) with the unit 
sphere C L* S , with respect to the metric g . 

Thus, in all cases S^ot turns out to be a manifold (see Section lL2l and |25]). □ 

4.6 Corollary. In the non degenerate case, the choice of a configuration rj-ot G Srot 
and of a scaled orthonormal basis in L* (g) S yield the diffeomorphisms (via the action of 
0{S,g) on Srot) 

S-rot — 0{S,g) ^ 0(3) , in the strongly non degenerate case; 
Srot — SO{S,g) ~ S0{3) , in the weakly non degenerate case. 

In the degenerate case, the continuous choice of an orientation on the straight lines 
(rrot) C S generated by each configuration r^^t ^ "S^rot and of a scaled orthonormal basis 
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in L* (8> S yields the diffeomorphisms 

Srot — S^(L* ® S,g) ~ S^(3) , in the degenerate case. □ 

From now on, in the non degenerate case, we shall refer only to one of the two con- 
nected components of Sj-ot , for the sake of simplicity and for physical reasons of continu- 
ity. Accordingly, we shall just refer to the non degenerate case (without specification of 
strongly or weakly non degenerate), or to the degenerate case. 

4.2 Tangent space of rotational space 

The rotational space S^ot is embedded into the environmental relative vector 
space Srei ■ Hence, the tangent vectors of Srot can be regarded as multi-vectors 
of Srei J which respect the rigid constraint. Here, we show a very geometric way 
to describe the tangent vectors of Srot , so obtaining a geometric interpretation of 
classical formulas of the mechanics of rigid systems. 

4.2.1 Non degenerate case 

First, we recall the following well-known result (see, for instance, |25j). 

4.7 Lemma. By regarding the orthogonal group as a subspace 0{S,g) C L{S, S) , 
we can identify the Lie algebra so{S, g) with the subspace so{S, g) C S* ® S , consisting 
of the tensors which are antisymmetric with respect to g . 

Proof. It is well-known that the subspace 0{S,g) C L{S,S) consists of the invertible elements 
/ € L{S, S) , which are invertible and fulfill the condition f^^ = f^- 

Hence, if c : 5? — *■ 0{S, g) C L{S, S) is a curve such that c(0) = id , then we obtain id ~ coc^^ = coc* , 
hence = D{c o c*)(0) = Dc{0) + Dc^{0) . Thus, the vectors tangent to 0{S, g) at the identity consist of 
antisymmetric endomorphisms. 

Conversely, we can prove that, for each antisymmetric endomorphism lo G L{S, S) , there is a curve 
c:R^ 0{S, g) C L{S, S) , such that c(0) = id and Dc(0) = lo . QED 

4.8 Proposition. We have the natural parallelising linear isomorphism 
(4.11) T,ot -T Srot S, ot X so{S,g) 

and the associated projection 

Prot : TSj-ot C Srot X S'rcl ^ So{S , g) : (rrot , ^^rot) ^ • 

The expression of the inverse isomorphism is 

''"rot • "^rot X S0{S,g) TS^ot C S^ot ^ S^d '■ {'l^rot, ^) ^ ^'"rot; 

Thus, for each (r^ot, "Wrot) ^ TS^ot C. S^ot x Sj-ei , there is a unique u G so{S,g) , such 
that Vi = Lj{ri) , for 1 < i < n . 



30 



4 Rigid systems 



Proof. By recalling that 5'rot is an affine space associated with the Lie group 0{S,g) , Lemma [T!^ 
yields the isomorphisms Trot ■ 

Next, in order to compute t^^I , let us consider an element uj G so{S,g) . Then, there exists a map 
a; : JR — !• SO{S, g) such that a)(0) = id e SO{S, g) and {Duj){0) = uj . Hence, for each rrot e 5rot C Srd , 
the curve c : Ft ^ 5'rci : A i-^ (a;(A)(ri), . . . ,a)(A)(r„)) is valued in Smt because 

\\Q{X){n)^Q{X){r,)\\ = MX)in - r,)\\ = \\n -TjW, VXeR, VI < ij <n. 

Hence, the tangent map Dc{0) £ is valued in Tr^^^S^ot G • On the other hand, we have 

-Dc(0) = (w(ri),...,w(r„)).QED 

Later, we shall give an explicit expression of the parallelisation Trot , via the "inertia 
isomorphism" ( Corollary 14 .241) . 

We can read the parallelisation Trot in a further interesting way, by means of an 
algebraic re-interpretation of so{S,g) . 

For this purpose, we recall the cross products X of S and of S* , defined by 

u X V := *{v A w) := g\i{u /\v)r]) = i(^g^{u) A g^{v)^ fj , \/ u,v ^ S , 
a X (3 := *{a A (3) := g\t{a A(3)f])= t[g\a) A g^{(3)) 7], Va, /? G 5* . 

The cross product commutes with the metric isomorphisms, i.e. we have 

g{u xv)= g\u) X g\v) and g\a X /?) = g^a) X gK(3) . 

For short, we set 

a Xmul f3mu\ := (a X . . . , a X /3„) , Va G 5* , V/3mui e -S'^ui • 
Moreover, we introduce the scaled vector space 

(4.12) Fang:=L-^®5. 

4.9 Lemma. The metric isomorphism (7^ : S*®S IJ^^S*^S* : a®v t— > a<^g{v, ■) 
and the Hodge isomorphism * : A'^S* h^'^ ^ S : 00 ^ i{uj) fj yield the linear 
isomorphism 

*o g^ : so{S, g) Vang : uj ^ i[g^{uj)^f] . □ 

4.10 Corollary. We have the natural parallelising isomorphism 

(4.13) Tang := * O g^ O T,ot ■■ TS, -ot ^ S^-ot X Vang 

and the associated projection 

Pang • TSj-Qt C Sj.Qt X Srei > Vang • ('"rot; '^^'rot) ^ ^ ■ 
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The expression of the inverse isomorphism t^^^ is 

^ang : Srot X Vang TSrot C 5rot X 5rel : {Vrot , ^) ^ (^rot , {^1 X n, . . . , Q X r„)) . 

Thus, for each {vj-ot, Vmt) £ TS^ot C S^ot x S^d , there is a unique Q E Vang , such that 
Vi = Q X Ti , for 1 < i < n . 

Proof. It follows from Proposition 14 . 81 and the Lemma lOl QED 

Later, we shall give an explicit expression of the parallelisation Tang , via the "inertia 
isomorphism" (Theorem I4.22|) . 

The above Corollary 14. 101 is just a geometric formulation of the well-known formula 
expressing the relative velocity of the particles of a rigid system through the angular 
velocity. 

4.11 Corollary. The transpose (T^^g)* of the isomorphism r^ng has the expression 
(4.14) (r-y* : T*5,ot ^ 5,ot x : (r,ot,«) ^ (r.ot, T.9\r,) X a,) . 

i 

Proof. The expression of t~^^ and cyclic permutations yield 

(^ang)*('-rot,«) {n) 1 = 

= a(Tang(»'rot,i^)) ^ a {fl X i„ui r,.ot) = X!"' ^ '^^^ X! (s'l"*) > X n)) 

i i 

i i i 

The cross product x is equivariant with respect to the action of SO{S,g) . Hence, 
the isomorphism Tang turns out to be equivariant with respect to this group. 

4.2.2 Degenerate case 

Let us consider the quotient vector bundle [Smt x so{S,g)] over S^ot of the vector 
bundle S^ot x so{S,g) , with respect to the vector subbundle /ifSmt] consisting, for each 
Trot € 'S'rot , of the isotropy Lie subalgebra /i[ri.ot] C so{S, g) of {r^-ot) C S (see, for instance. 

Now, let us refer to the degenerate case. We can rephrase the results concerning the 
non degenerate case by a quotient procedure. In particular, we have the following results. 

For each £ S^ot , the isotropy Lie subalgebra associated with Trot consists of 
antisymmetric endomorphisms (f) E S* ^ S which preserve the 1-dimensional vector 
subspace (rj-ot) C S generated by Trot • 
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4.12 Lemma. We have the natural hnear fibred isomorphism 

] : T5rot [Srot X so{S,g)] . 

Proof. It follows from a well-known result on homogeneous spaces (see, for instance, |^). QED 

Let us consider the quotient vector bundle [Srot x ^ang] over Srot of the vector bundle 
Srot X ^ang with respect to the vector subbundle a [Srot] consisting, for each Trot G Srot , 
by the 1-dimensional vector subspace (rrot) C Vang generated by Trot ■ 

4.13 Proposition. We have the linear fibred isomorphism 

''"ang] • ^Srot ^ [Srot X V^ang] • 

The expression of the inverse isomorphism [r^^g] is 

''"ang] • [Srot X V^ang] ^ ^Srot C Srot X S^ul '■ 

■ (rrot , [^rot, ^]) ^ [rj-ot , X ri, . . . ,Q X r„)) , 

where the cross products Qxri turn out to be independent of the choice of representative 
for the class [rrot, ^] • 

Thus, for each (rrot,^^rot) ^ TSj-ot C Srot x Srei , there is a unique [rrot,^] G [Srot x 

Vang] rrot 5 SUch that Vi = Q X Vi , ioT 1 < i < u . U 

4.14 Proposition. A continuous choice of an orientation of the straight lines (rrot) C 
S generated by the configurations rrot € Srot yields the linear isomorphism 

TSrot ^ TS\h* S,g).n 
4.3 Rigid system metrics 

The multi-dynamical metric of Smui induces a metric on Srot , which can be 
regarded also in another useful way through the isomorphism Tang , and will be 
interpreted as the inertia tensor. 

Moreover, the standard pattern metric of Vang induces a further metric on Srot • 

4.15 Proposition. The inclusion i^g : Prig ^ -Pmui yields the geometrical and 
weighted scaled Riemannian metrics 

(y'rig •= Vig fi'mul • -F*rig IL ® (T -Prig ® T Prig) ; 
Grig •= "^rig ^mul • -frig ^ L. CS) (T Prig ® T Prig) • 

The splitting Prig = Pdia © Srot is orthogonal with respect to the metric Grig • 
Proof. The splitting Pmui — Pdm © S'l-ci is orthogonal with respect to the metric Gmui and we have 

-Prig = Pdia ® Srot , with Srot C S re\ ■ QED 
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4.16 Proposition. The inclusion i^ot '■ -S^rot ^ •S'rei yields the geometrical and 
weighted scaled Riemannian metrics 

fl'rot ■— Cotfl'rel • ^rot — ^ C?) (T*/S'rot 2^*'S'rot) ) 
Crot ■= Cot ^rel : S'rot ^ ® {T*S rot T*S rot) ■ 

For each (fi Xmui ^rot) , Xmui ^rot) ^ ^rrot-^rot C S^ei , we have the expressions 

5'rot(?^rot) Xmui rrot , ^' X^ul ^rot) = 1] ^(^' ~ ^(^'' ^) ^')) > 

i 

Grot('^rot) Xmui rrot , ^' Xmui ^rot) = 1] /^i (s'l'^i, '^i) fl'(^, " 5'('^i, ^) 5'('^i, ^O) • 

i 

Proof. In virtue of standard properties of the cross product, we obtain 

ffrot(''rot) ^mu\ ''rot , ^' ^ mu\ ''rot) — ^9{^ X Tj, fl' X Tj) =^g[{n' X n) X Q , Tj) 

i i 

= ^ (fl(ri, Vi) g{n, fl') - gin, Q) g{n, n')) . QED 



We can regard the metrics g^ot and Grot in another interesting way, via T; 



ang ■ 



4.17 CoroUciry. In the non degenerate case, the isomorphism Tang allows us to read 
g^rot and Grot as the scaled metrics 

(4.15) a := (r-„y * ^^rot : ^^ot ^ ® {Vl^^ ® VIJ , 

(4.16) S := (r-,y * Grot : 5rot - ® {Vl,^ ® , 

with expressions 

a{rrotm, = E {air,^ ^0 ^(^, - ^(^., ^) ^(^^, ^O) , 
E(rrot)(^^, ^^') = E /^^ (^(^i' ^(^' - 9(ri, n) giu, Q')) . □ 

i 

We have a further natural metric of Srot . For this purpose, let us consider the metric 
g G V*jjg (g) V*^g of Vang naturally induced by the pattern metric gi of . We can make 
the natural identification 0{Vi^^g,g) ~ 0{S,g) . 

4.18 Proposition. In the non degenerate case, we obtain the unsealed Riemannian 
metric 

fl'ang •= Tang 9 • "^rot ~^ T*Srot ® T*Srot ■ 
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For each {r^^t , ^ Xmni ^^rot) , ('"rot , ^' Xmul Trot) € TS j-ot C S^ot X 5rei , we have the 
expression 

5'ang('"rot) Xmui Trot , ^ Xmul ^rot) = • 1^ 

All metrics of Srot considered above are invariant with respect to the left action of 
0{S,g). 

4.19 Proposition. In the non degenerate case, the choice of a configuration Trot G 
Sj-ot and of an orthonormal basis in Vang , respectively, yields the following diffeomor- 
phisms (via the action of SO{Vi^^g, g) on Sj-ot) 

Srot ^ SO{V,^^,g) ~ ^0(3), 

which turn out to be isometries with respect to the Riemannian metrics gang, — |^ang 
and —^k^, of Srot , ^ang and SO (3) , where k^ng and k^ are the Killing forms. 

Proof. The above diffeomorphisms yield the Unear fibred isomorphisms 

Tr„tSrot ^ S0{Ve,ng,g) ^ So(3) . 

On the other hand, the natural isomorphism so(Vang,5) ^ ^ang , induced by and ^ , is metric. 
Hence, in virtue of the definition of (jang , the isomorphism Tr^^^^S'rot — so(Vang, 5) turns out to be metric. 

Moreover, the metric g of Vang turns out to coincide with the metric — ^ fcang of so(V^ang, 5) • In fact, 
we have g{uj,uj') = —\ tr ((ilx) o (fi'x)) . It is easy to see that the isomorphism so{Va.ag,g) — so(3) is 
metric. QED 

We leave to the reader the task to extend the above results to the degenerate case. 
4.3.1 Angular automorphisms 

4.20 Proposition. The unsealed metric g of Vang allows us to regard the metrics a 
and S as scaled symmetric fibred endomorphisms 

a:=g^oa' : S^ot ^ ® (K.^ ® Fang) , 
t:=g^ot' : S^ot ^ ® {Vl,^ ® Fang) • 

We have the expressions 

i i 

S(rrot) (fi) = ^^iiri {VL Ti) fii [g{ri, r^) Vt - g{ri, Q) r^) . 

i i 

In the non degenerate case, the above maps are automorphisms. 
Proof. It follows immediately from CoroUarv l4.17l QED 
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4.21 Lemma. In the non degenerate case, we have 

= cr" o t?^ and E"^ = E« o c?^ 
(a-i)* = ^^oa« and (S-^)* o E« . 

Proof. Wc have (j-i := (g« o cr'')-i = (cr^)-^ o (^tt)"! = cr" o c,^ 

Moreover, the symmetry of g and a give (a"^)* = (tr" o gi'')* = (5^)* o (cr")* = 51'' o a" . 

Analogous proof holds for S . QED 

The automorphisms a and S yield the following explicit expressions of the map Tang . 

4.22 Theorem. In the non degenerate case, the isomorphism Tang has the expression 

"Tang ■ TSrot C S-rot ^ 'S'rel ^ 'S'rot ^ Vang : (r-rot , ^^rot) ^ (^rot , ^) , 

where 

(4.17) n = a-\r,ot){Yl^i X ""0 = ^-^(rrot)(l]Afi n X . 

i i 
Proof. Let rrot e Srot , ^rot = {vi,. ..,Vn)e Tr^tS'iol C Srel and set := Pang(?'rot, l^rot) e Vang • 

The definitions of a and of grot yield, respectively, the following equalities, for each e Vang , 

5rot(rrot) (t^rot , f^' X mul ''rot) — (^irrot) , 

5rot(rrot) {Vrot , ^' X mul ''rot) — ^ 9{Vi , O.' X Tj) = g{Q' , X Vi) = C/(^ X t^j , O') . 

Then, by comparison of the above equalities, wc obtain cr''rot(^'rot)(f^) = .9^(X]i^« X Vi) , hence 
CT(rrot)(fi) := (fl" o a^iot){riot){^) = Y,in X Ui , which yields O = c^^^{r^ot){Y.t X 1^2) • 

We can prove the 2nd expression of Vt in analogous way, by replacing grot with Grot • QED 

In the classical literature, Q is computed by means of the Poisson's formulas, in terms 
of a basis. The above Theorem provides an intrinsic expression of Vt , which plays an 
essential role in next sections. 

4.23 Note. The map Tang is a geometric object, which has nothing to do with masses 
and weights, because the rigid constraint does not involve the masses. Accordingly, the 
1st formula in the above Theorem is natural, while the 2nd one sounds quite strange. 
Indeed, the 2nd formula is true for any arbitrary choice of the weights. We have added 
the 2nd formula for the sake of completeness. Actually, in the 2nd formula, the weights 
appear both in the expressions of the sum and of E ; eventually, the contribution of the 
weights disappear. 

In order to help understanding this situation, we prove directly that the 1st formula 
implies the 2nd one. 

Let Trot £ "S^rot ) '^^rot ^ ^^rrot'^rot C S^^l ^-^d Set Q := Pang(^rot ) "^^rot) ^ ^ang • 
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Then, the definition of S and the 1st expression of pang imply 

t{r,ot)i^) = ^/iiTi X (^] X ri) = Y,l^iriX (('^r"ot(^rot)(5I X ^i)) X ^i) • 

i i i 

On the other hand, in virtue of the definition of pang and of its 1st expression, we have 

(^mi('^rot)(^ri X Vi)) Xri = Vi. 
i 

Hence, we obtain S(rrot)(^) = J2i A^i ^iXVi, which yields Q = ^~^{rj.ot){J2i f^i fiX'^i) ■ □ 

4.24 Corollary. In the non degenerate case, the isomorphism Trot has the equivalent 
expression 

Trot : TSrot C Srot X Srcl ^ Srot X (5* S) : {Vrot , t^rot) ^ (rj-ot , , 

where 

i i 
Proof. It follows from Proposition 14. 221 and from the composition of algebraic isomorphisms 

y^„g = L-i ® S L2 ® A^S* S*(^S. QED 

4.25 Proposition. The eigenvalues of S turn out to be constant with respect to 
Srot , in virtue of the invariance of S with respect to SO{S, g) . 

In the non degenerate case, we have three eigenvalues. Then, three cases may occur: 

A := Ai = A2 = A3 , spherical case, 
A := Ai = A2 7^ A3 , symmetric case, 
Ai 7^ A2 7^ A3 7^ Ai , asymmetric case. 

In the degenerate case, we have two coinciding eigenvalues 

A := Ai = A2 = ^pj5'(ri,ri) . 

i 

Analogous results hold for a . 

We have studied the diagonalisation of S with respect io g . In an analogous way, 
we can diagonalise Grot with respect to Gang • Indeed, in this way we obtain the same 
eigenvalues and the same classification, because the two diagonalisations are related by 
the isomorphism Tang • 

The scaled metric mo S , or the scaled automorphism tuq E , are called the inertia 
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tensor and the scaled eigenvalues Jj = mo Aj : Sj-ot — >■ (L^ <8) M) ® iR of the inertia tensor 
are called principal inertia momenta. 

4.3.2 Continuous interpretation 

We can interpret the above results concerning the parallelisation of /Srot also 
in terms of continuous transformations. Here, in order to keep the thread of our 

reasoning, we adopt a purely geometric approach which does not involve time, but 
this section can be easily rephrased in a true kinematical way, by replacing IR with 
T , or T (g) iR , as appropriate. 

We define a continuous transformation as a map 

C:]Rx{RxP)^P, 

such that, for each t,t' ,t ^ IR, p e P , 

C{0,t,p)=p and C{T',t + r,C{r,t,p)) ^C(r + T',t,p). 

A continuous transformation is said to be rigid if, for each T,t & M, p,q & P , 

\\Cir,t,q)-Cir,t,p)\\^\\q-p\\. 

We can prove that a continuous transformation C is rigid if and only if its expression 
is of the type 

C(t, t,p) = c{t) + ip-o), yT,teM,peP, 

where a e P , c : M ^ P and ^ : M x M ^ SO{S, g) . 

Let us suppose that C be rigid. The partial derivative of $ with respect to time, at 
T = , turns out to be an antisymmetric endomorphism 

5^ ■.lR^so{S,g) C 5* 

Hence, the velocity of the continuous transformation v : T x P — > is given by 

v{t,p) = Dc{t) + 6^t) (p-o), yt eT , p e P . 

On the other hand, we obtain the map 

Q :=(*oc;^)((5$) : IR ^ L* ® S , 

via the metric isomorphism g'" : S S* and the Hodge isomorphism * : A'^S* ^ h* S . 

Therefore, we can express the velocity of the continuous transformation by the classical 
formula 

v{t,p) = Dc{t) + n{t) x{p-o), yte pe P. 
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4.26 Note. Let P^g C Pi x . . . P„ be a non degenerate rigid configuration space 
and Srig : IR — > Prig be a map. 

Then, there is a unique continuous rigid transformation such that the particles of the 
continuous transformation, which coincide with the particles of the discrete rigid system 
at a certain time, move as the particles of the discrete rigid system. 

In other words, there is a unique rigid continuous transformation 

C : iR X (iR X P) ^ P , 
such that, for each p = (pi, . . . e Prfg , 

Then, for each p — (pi, . . . G P^g and t & M, we have 

V{t,pi) = dSi{t) = dScen{t) + ^{t) X . 

Indeed, the rotational components of the velocity of the continuous and discrete rigid 
maps coincide. □ 

4.4 Splitting of the tangent and cotangent multi— space 

Next, we study the relation between the tangent and cotangent spaces of the 
rigid configuration space and the tangent and cotangent spaces of the environmen- 
tal space. Namely, we exhibit a natural orthogonal splitting of the environmental 

tangent and cotangent spaces into three components: the component of the center 
of mass, the angular component and the component orthogonal to the rigid con- 
figuration space. This splitting will have a fundamental role in mechanics of rigid 
systems. 

4.4.1 Splitting of the tangent multi— space 

Let us consider the space TPmuilp^.^ = Prig x Sj-^i ■ 

4.27 Theorem. We have the orthogonal splitting, with respect to Gmui 

^-Pmul|Prig ~ TPj-ig © T"'"Prig = (TPcen X TSj-qi) © T^P^{g , 

Prig -Prig 

where T^P^g is the orthogonal complement ofTPj-ig in TPmui|p^;g . 
The suhspace T-^P^g is characterised by the following equality 

T^Prig = {(Prig , ^^mul) £ -Prig X 5niul | X] A*^ = ^ ' Z] '^^ ^» = 0} ^-PmullPrfg ■ 

i i 
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Moreover, the expressions of the projections associated with the splitting are 



(4.18) 


"cen 


^-Pmul Prig " 


ccn 


• (Prig ) 


^^mul) ^ 


(Pcen ; 


^ccn) 5 


(4.19) 




^-Pmul Prig " 


TSrot 


• (Prig ) 


^^mul) ^- 


{Trot , 




(4.20) 


vr^ 

rig 


T'-Pmul Prig ~ 


-f ± -IT rig 


• (Prig ; 


^'mul) ^ 


(Prig , 


^mul ^dia X j^mj ^rot ) ; 


where 

















Pcen = O + '^HiiPi - a) , Pdia — (Pcen, • • • , Pcen) , T^ot = Prig - Pdi 

i 

Veen ^ ^ A'"* ) ^dia • (^cen; ■ ■ ■ ; ^cen) ; ^rot ^mul ^di 

i 



la ; 



aia ; 



Proof. The expression of the 1st projection is obvious. 

Let us prove the expression of the 2nd projection. For each r^ot G •S'rot i ''^rci G S,-e\ and fl' S Vang j 
in virtue of the definitions of S and of Gici , and by a cychc permutation, we obtain the equahties 

Grcl {Vicl , f^' Xmul '"rot) = Grol (TvTrot (^ot , f rcl) , ^' X mul ^ot) = Grcl X rrm\ fiot , ^' X- rrm\ ''rot) 

Grcl (I'rol , ^' Xmul ''rot) = ^ Sl^^i - ^' X = ^(ri' , ^ /i^ X 1;^) = .g(^ /i^ n X Vi , fl') . 

i i i 

A comparison of the above equahties yields E''(ri-ot) {^) — ^''(Si /^i ''» X Wi) , hence 
fl = (S''(rrot) ° g^)C^^J■iri X Vi) = S"^(ri.ot)(^M»'"j X Vi) . 

i i 

Then, the characterization of T^Prig is easily obtained by considering the multivectors whose pre- 
vious projections vanish and by recalling that Tang and S are isomorphisms. 

Eventually, the 3rd projection is obtained by subtracting the previous projections. QED 

We observe that the expression of TTTrot is similar to the 2nd formula of Proposition 
14.221 However, we stress that the multivector Wmui in the above Theorem needs not to 
be tangent to S^t and its projection on S^ot involves the weights. In the particular case 
when the multivector v^ani is tangent to S^t , the expression of TTTrot reduces to the 2nd 
formula of Proposition 14.221 

4.4.2 Splitting of the cotangent space 

Let us consider the space T*Pmui|Prig •= -Prig x S*^^ . 

4.28 Theorem. We have the orthogonal splitting, with respect to Gmui 

rp* ID I n^* C! m T'* p ('T'* / J v T'* C \ a\ T'* p 

-f^ mul I Prig — '-'rig ty -f_|_J^rig — yJ- -f^ccn ^ »-'rotJ W -t_|_-rrig; 

Prie Prie 
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where T^Pj-ig is the orthogonal complement ofT*Pj.[g in T*Pniui|Prig (i-^-, the space of 
annihilators ofTPj-ig). 

The subspace TJ^Prig is characterised by the following equality 



rp* rf 

-L rig 



Prig ) Oi^xil, 



e Prig X I ^fti = , Y.9\vi) X = O} C T*P^,,i\p^^^ . 



Moreover, the expressions of the projections associated with the splitting are 



(4.21) 
(4.22) 
(4.23) 

where 



" cen • ^ mul | Prig con 

T* rp* I y I rp* Q 

TTrot • J mul I Prig — ^ ^ '-'rot 

7r_L : T*Pjnul|pHg ~^ rig 



(Prig ) '-'ijjiul) ' (Pcen j ^cea) 

(Prig , ftmul) ^-> ('^rot , a^ot) 

iPng , ttmul) ^ {rrot , "mul " "dia " "rot) 



Pc 



(4.24) 



i 

i 



Pdia • (Pcen; • • • y Pcen) ; '"rot Prig Pdia ; 

Q^dia (/^l C^cen; • • • ; f^n C^cen) ; 

^{J29\ri) X OiU XmuiG^uil'-rot). 



Proof. The commutative diagram 



..4 



mul I Prig 



T*irro 



Theorem 14 . 2 71 and the definition of E give 



) = (GLoT7r„toGLl)(Prig,amul) 

= (Grot °T'7^rot)(prig, ^17 ff" ("l )'■■•: 7^ (««)) 



mul ''"rot 



= (w, ((<?^oS«)(r,ot)(E5'('^0 X a,)) X„,,aGUi(r-rot)^ 

i 

= (r,.ot, ((S^i)*(r,ot)(^5'(?-.) X a,)) X„,,i G^,,i(r,ot)) . QED 



The projection T*7ri.ot can be expressed in terms of Vang • In this way, we recover 
the classical formula of the "total momentum" of a multi-form. Here, this formula arises 
naturally from our geometric interpretation of T*S'rot • 
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4.29 Corollary. We have the projection 



(4.25) := (r-„y * o TV^ot : T*P^,i|p,^ ^ S^ot x V^,^ 



■ (Prig, ttmul) (^rot, XlS'^ln) X 

i 

where r^t := TTrotbrig) • 

Proof. By recalling the expressions of T*7rrot , , (Ta^g)* and S , we obtain 

<Sang(Prig, Q!mul) = ((Ta^g)* O TVrot) (Prig, amul) 

= (ra-ng)*(r-rot, ((S"^)* ) (^i ) >< X G^.^! )) 

= (r.ot, X (((S-')*(rrot)(E5'(^i) X a,)) X Mifl'(n))) 

i 3 

= (vrot, (a^^^'i) Xmul G^l (S"' (rrot) (J^ X ^^K))) X r,ot) ,)) 

i 3 

= (r-rot, g^^iJ-in X (s-^(rrot)(X] ''j' X aKoii)) x n))) 

= (j-rot, ^<?^(n) XQi).QED 



4.5 Kinetic energy and momentum of the rigid system 

According to our scheme, we define the rigid kinetic energy, the rigid kinetic momen- 
tum, the rotational kinetic energy and the rotational kinetic momentum as 



A^rig 


— ^rig^mul 


T" 




5 TPrig - 


(T-2 




5M) C 




^^rig 


imoGrigKig, 


f^rig) ) 


'Prig 


^rig^mul 


T" 




5 TPrig - 




(g) 






t'rig 1- 


mo Grig(^rig) , 




^rot 


*rot'^mul 


T" 




5 T5rot - 












|moGrot(^'rot 






•~ ^rot^mul 


T- 




5 T5rot - 




(8)L^ $ 


§M) $ 






-^^0 G'^ot(^rot)- 





Then, by taking into account the angular parallehsation Tang : TS^-ot —>■ -S^rot x V 



ang ) 



we obtain the angular kinetic energy and the angular kinetic momentum 



ang 



T*ng ^rot : ® TS,ot ^ (T'^ (g) M) ® iR 



ang 



rot 



-1 



V 



ang 



: (?^rot, t^rot) ^ mo E^(rrot)(f^) = "^0 X!'"^ (din, Ti) g\0) - g{ri, g^ 
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where Q := Pang(^^rot) ■ 

According to the sphttings TP^g = TPcen x TS^t and TPj-ig = TP cen x (S'rot x ^ang) , 
we have, respectively, the sphttings 

^rig ^con ~l~ ^rot ; "^rig ("^ceri) T^rot) ; 

^rig ^ccn ~l~ ^ang ; ^rig (^ccrij ^ang) • 



4.6 Connections induced on the rigid system 

First of all, in the non degenerate case, the generalised affine structure of P^g induces 
a fiat connection Vafr (see Section 

Moreover, according to the Gauss' Theorem (see, for instance, (Tj), the geometric 
metric g^i^ and the weighted metric Grig induce two distinct connections Vrig and VRig , 
respectively, on P^g . Actually, we shall be mainly concerned with VRig , which is the 
most important of the two, because of its role in dynamics. 

4.30 Proposition. The connection VRig splits into the cartesian product of the con- 
nections Veen and Vrot of Peen and St-oi . □ 

4.7 Kinematics of rigid systems 

In this section, we apply the splitting of TPmuilPHg to the velocity and the 
acceleration of a rigid system. Namely, the velocity splits into the two components 
of the center of mass and the velocity relative to the center of mass. On the other 
hand, the acceleration splits into the three components of the center of mass, relative 
to the center of mass and the term given by the 2nd fundamental form of the rigid 
configuration space. 

Let us consider a rigid motion s^g '■ T P^g . 

4.31 Proposition. According to Corollary I4.1UI and Theorem 14.221 we obtain the 
splittings 

"^rig ("^cenj '^[■ot) • T >■ Peen X Sj-ot 

fi)) : T ^ (T ^ (g) TPccn) X (^Srot X (T ^ (g) Vang)) , 

where 

Seen ■— "^ccn ° -^rig ; ■^j-ot •= ''Trot ° ^rig ; ^2 .= S (^^ /ij (Si.ot)j X d(^Sj.Qt)i) . □ 

i 

The map Q := Pang ° ^^ang ° c^-^rig : T ^ T* eg) Vang is called the angular velocity of 
the rigid motion. 
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4.32 Theorem. The acceleration splits into the three components as 

+ VrotC^Srot + N O dSrot , 

where N : TPrig —>■ T-^Pj-ig is the 2nd fundamental form of the connection Vmui , with 
respect to the metric Gmui • We have the expressions 

(4.26) Pang(VrotC?Srot) = d^l + ±-\s,ot) X ±{Srot){^)) , 

(4.27) N{s,ig){Q) = Qx^^i{Qx 

mul "^rot) ^ ("5 rot mul -^rot • 

Proof. The proof is analogous to that for the case of one constrained particle. 

In virtue of Theorem 14. 2 21 Corollary 14 . 1 Ul and the Leibnitz identity for the cross product, we obtain 

Pang(Vrot(^Srot) = 

= (Pang O r7rrot)(VmulC^Srot) 

= S"^(Si.ot)(^Ali(Srot)i X Vd(Srot)i) 

i 

= S"^(Srot)(^/^i(Srot)i X V(fi X (Srot)j)) 

i 

= S"Hsrot)(Xl/^»(*™t)» ^ ('^^ ^ (Srot)*)) + (Srot ) ( ^1 ('5™* X (f7 X (SI X (s,.ot),))) 

i i 

= t-^{Srot){^f^z{Srot)t X {dil X (Srot)^)) + {Srot) f^t^ X {{s,ot)t X {ft X {s,ot)z)) 

i i 

= dn + ±-\s,ot){n X E(s,ot)(r!)) . 

Moreover, we have 

N{dSrot) = 

^mul'^^rot ^rot^-^rot 

= dD, Xjuui Sj-ot + ft X,uui {n X,„ui Slot) — Xmui Srot ~ S ^(Srot)(S^ X S (Sfot ) (S^)) X ,„ui Sjot 
= nX {n X„iui Srot) - S"^(Srot)(f^ X E(Srot)(f^)) X S^ot • QED 

4.33 Note. The map dfl is the covariant derivative of Q with respect to the natural 
connection Vas of Srot induced by Tang - Hence, the map E~^(srot)(^ X S(si.ot)(^)) ex- 
presses the Christoffel symbol of the connection Vrot with respect to the parallelisation 

''"ang ■ n 

4.8 Dynamics of rigid systems 

In this section we study the equation of motion for a rigid system. 

According to the results of the previous section, we show that the equation of 
motion in the environmental space splits into three components: the equation of 
motion for the center of mass (related to the linear momentum), the equation of 
motion for the relative motion (related to the angular momentum) and equation for 
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the reaction force (related to the 2nd fundamental form of the rigid configuration 
space). 

4.8.1 Splitting of multi-forces 

According to the scheme discussed for a constrained system, we assume a multi-force 



^mul : JlPmnl ^ (T"^ ® ® M) ® T* P 



mul 

and consider its restriction to the phase space of the rigid system 

F := F^ui I JiPh, : -^i-Prig ^ (T'^ ® ® M) ® T*P^,i . 

4.34 Proposition. According to the splitting of T*i-'mui|-Prig ? we can write 

-1 -1 rig -1 rig ± ; 

where 

^rig = ^rig^ ■ JiPrig (T^^ ® ® M) O T*P . 
Moreover, according to the splitting of T*P.^i^ (see Theorem I4.28|) . we can write 

-^rig (-^ccn; -^rot) (-^cenj -^ang) ) 

where F^en and Fang turn out to be, respectively, the total force and the total momentum 
of the force. We have the following expressions 

F =T*-n oF- =^F- F =S oF- =^r xF 

^ cen "cen ^ ^ rig / , ^ t i ^ ang '-'ang ^ ^ rig / , Li j. j , 

i i 

Frot ((^ ) (-^ang)) X jjiul Z^rig ; -^dia (A''! -^cen; • • • ; /^n -^cen) ; 

FrcoA -^dia -^rot ; 

where we have set Trig : 5rot C 5roi ^ (g) S*^^ : Trot ^ G'^^i(rrot) • □ 

4.35 Note. If the multi-force Fmui fulfills the 3rd Newton's principle then its compo- 
nent tangent to the constraint Frig vanishes. 

Proof. If Fmui fulfills the 3rd Newton's principle, then we obtain i^con = and -Fang = 0, hence 
Frig = . QED 

Moreover, we assume a reaction force 

R : JiPrig ^ ® (g) M (g) T*P^^x , 

which splits analogously to the force. 

The above Note holds also for the reaction. So, we assume i?rig = , i.e. R = Rj_ . 
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4.8.2 Splitting of the equation of motion 

Eventually, we are ready to split the equation of motion into three components. We 
follow the scheme developed for one constrained particle with the additional results arising 
from the present framework. 

4.36 Corollary. The rigid motion s^g : T Pj-^g and the reaction force R fulfill the 

Newton's law of motion 

if and only if 

(4.28) mo Gcen(VcenC?Scen) = -^cen O jiSng 

(4.29) mo Gl^,{N o ds.ot) 
where 

Pang(VrotC?Si.ot) =(1^1 + S"^(Srot)(fi X S(Srot)(f^)) • 

Proof. The proof follows from Theorem QED 

The angular component of the above equation of motion is referred to as Euler's 
equation. 

4.37 Note. Contrary to the general constrained case, here the equation on the con- 
strained space can be split into the center of mass and rotational components, due to the 
fact that the rigid constraint is without interference between particles and the center of 
mass. 

On the other hand, as in the case of a system of free particles, we cannot solve the 
first two equations independently, unless the total force factors through the projection 

TTcen OU P cen ■ □ 

4.38 Corollary. The reaction R± is given by 

(4.30) = X (O X Trig) - S-l*(Srot) ((O X S^(s,ot)(l^)) X^,i Trig) 

-^mul ~l~ (pi Pceni • • • i Pn -^cen) ~l~ ((^ ) ("^a^ng)) Xmul ?^rig • 

where we set := g^{VL) . 

Proof. The reaction R± is determined by on solutions of the equations of motion by the following 



^0 ^^Tang(Vrot'^'5rot)) — -^ang O Jl^rig 

(Fj_ + R^) o jiSrig , 
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equalities 

-R_L = G^ul(^ ° dSrot) - F±0 jiSng 

= G^ui(Ar o dSrot) - [Fmul - -Fdia " F^ot) O JlSrig 

mul ^rot) 

JlSrig Xmul ^rig 

Xmul Srig) 

Srig -Pcenj • • • ) A*n -^cen) ° Jl^rig jl ^rig X mul ^rig • QED 



Now, we express the Newton's law in Lagrangian form, in our special case of rigid 
systems. To this aim, we introduce an appropriate chart on Prig • We refer to a chart (x*) 
on Pcen and to a chart (for instance, the Euler's angles) ia^) on S'rot • Then, the induced 
chart on TQ is (x*, , x\ a^) . 

Suppose that Frig : -Prig (T^^ (g)L^ (8)M) (8)T*Prig be a conservative positional force, 
with potential Ung : Prig ^ (T^^ (g) (g) M) (g) I? , and that ii!rig = . Then, the induced 
Lagrangian function turns out to be the map 

Aig /Crig - Wrig : TPrig ^ ® ® M) ® iR . 

4.39 Corollary. Let Srig • T Prig be a motion. Then, Srig and the reaction force 
R fulfill the Newton's law of motion if and only if the following equations hold 

D [-of O ^-ng j - ^ ° ds,., = , D^-^O ds.,, j-^O Ci.rig = , 

GluliN O dSrig) = (Fi + R^) O jiSrig . □ 
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